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Abstract 

In 1,2 or 3 dimensions a scalar wave excited by a non-negative source 
in a viscoelastic medium with a non-negative relaxation spectrum or a 
Newtonian response or both combined inherits the sign of the source. 
The key assumption is a constitutive relation which involves the sum 
of a Newtonian viscosity term and a memory term with a completely 
monotone relaxation kernel. In higher-dimensional spaces this result holds 
for sufficiently regular sources. Two positivity results for vector-valued 
wave fields including isotropic viscoelasticity are also obtained. Positivity 
is also shown to hold under weakened hypotheses. 

Keywords: viscoelasticity, Newtonian viscosity, acoustics, relaxation, com- 
pletely monotone function, complete Bernstein function, MacDonald function 
Notation. 

[a, b [ := {x 6 K | a < x < b}\ 
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I: unit matrix; 
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£ - the cone of LICM functions; 



1 



25 - the cone of Bernstein functions; 

5 - the cone of complete Bernstein functions. 

1 Introduction 

Acoustic signals in real elastic and viscoelastic media reproduce the shape of 
the source pulse. In particular, a non-negative source pulse gives rise to a 
non-negative signal. This property is important in acoustic and seismological 
applications. We shall show that a large class of mathematical models of vis- 
coelastic media, including Newtonian viscosity and hereditary viscosity with 
positive relaxation spectrum, have this property. A viscoelastic model obtained 
by connecting elastic, Newtonian and Debye elements in parallel also has this 
property. Positivity of viscoelastic signals can be considered as a confirmation 
of the presence of these components. Unfortunately, it follows from a remark 
in Sec. [TU] that a more general class of equations have non-negative Green's 
functions and some of them might be compatible with viscoelasticity. 

Positivity of viscoelastic pulses was studied in a paper by Duff [8]. Duff 
assumed various special models with rational complex moduli. Duff's models 
are however loosely related to viscoelasticity and his assumptions are either 
excessively restrictive or inconsistent with viscoelasticity. 

In Sees 2HS] a general scalar viscoelastic medium with the constitutive equa- 
tion a = a e + G(t) * e with a completely monotone relaxation modulus G and a 
non-negative Newtonian viscosity coefficient a is studied. We show that a scalar 
viscoelastic wave field propagating in a d-dimcnsional medium and excited by 
a non- negative pulse is also non-negative provided d < 3. For higher dimen- 
sions and for non-zero initial data only wave fields excited by sufficiently regular 
sources are non-negative. 

Positivity can be extended to matrix-valued fields, e.g. to Green's functions 
of systems of PDEs. In Sec. [6] we consider a system of PDEs resembling the 
equations of motion of viscoelasticity with a CM relaxation kernel and prove that 
the Green's function of this system of equations is positive semi-definite. This 
result does not apply to general viscoelastic Green's functions, which involve 
double gradients of positive semi-definite functions. The signals propagating in 
isotropic viscoelastic media can however be expressed in terms of scalar fields 
representing the amplitudes or the potentials. It is thus again possible to prove 
that the scalar amplitudes excited by a non-negative force field are non-negative 
(Sec.©. 

In Sec. 19.11 positivity is proved under somewhat weaker assumptions, in- 
consistent with viscoelastic constitutive equations. In Sec. 19.21 we examine an 
example of a scalar differential equation which has a positive Green's function 
even when the assumptions of viscoelasticity with a positive relaxation spectrum 
fail to be satisfied. 

Our method takes advantage of certain classes of functions defined in terms 
of integral transforms and, alternatively, in terms of analytic properties of their 
complex analyitic continuations. Useful relations between these classes can be 
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established using the integral transforms. Analytic properties yield very useful 
nonlinear relations between these classes. These classes of functions appear in 
the viscoelastic constitutive equations [26] but they also turn out very useful in 
studying positivity of Green's functions. 

2 Statement of the problem 

In a hereditary or Newtonian linear viscoelastic medium a scalar field excited 
by positive source is non-negative. This applies to displacements in pure shear 
or to scalar displacement potentials. The key assumption about the material 
properties of the medium is a positive relaxation spectrum. The result holds for 
arbitrary spatial dimension. 
We consider the problem: 

P D 2 u = aV 2 Du + G(t)*V 2 Du + s(t,x) t>0, xeR d (1) 

with s(t, x) — 9{t) (cj + c 2 t) 8{x) and the initial condition 

u(0, x) = uo S(x), Du(0, x) = iio S(x), (2) 

(Problem I) as well as s(t, x) — c6(t) S(x) with a solution assumed to vanish for 
t < (Problem II). It is assumed that a > and G is a completely monotone 
(CM) function. 

The Laplace transform 

/•OO 

u(p,x):= e- pt u{t,x)dt, Rep>0, x G R d (3) 
Jo 

satisfies the equation 

pp 2 u(p, x) = Q(p) V 2 u(p, x) + g(p) S(x) (4) 

where 

Q(p) ;=ap + pG{p) (5) 

In viscoelastic problems a is the Newtonian viscosity. The instantaneous 
elastic response is represented by Go := lim p ^ 00 [pG(p)]. The equilibrium re- 
sponse is represented by := lim p ^ a [p G(p)]. Purely elastic response corre- 
sponds to Q(p) — const. 

The function g is defined by the equation 

9i(p) = ~ +pu n + u (6) 

in Problem I and 

9i lip) = 1 (7) 

in Problem II. 
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3 Basic mathematical tools. 



The classes of functions appropriate for viscoelastic responses are reviewed in 
detail by Seredyhska and Hanyga (author?) [26|, 15|. 

Theorem 3.1 

If the function u(-, x) is completely monotone for every x G K d , then u(t, x) > 
for every t > and x € R d . 

Proof. If u(-,x) is completely monotone, then, in view of Bernstein's theorem 
(Theorem l3.2l) . for every x e R d it is the Laplace transform of a positive Radon 
measure m x : 

u(-,x) = f e- ps m x (ds) (8) 

J[0,oo[ 

The Radon measure m x is uniquely determined by u(-,x), hence m x (dt) = 
u(t,x) dt is a positive Radon measure. Hence, in view of continuity of u(-,x), 
we have the inequality u(t, x) > for aU t > and x 6 M. d . 

□ 



The problem of proving that u(t, x) is non-negative is thus reduced to proving 
that u(-,x) is completely monotone. The crucial step here is the realization 
that Q in © is a complete Bernstein function. We shall therefore recall some 
facts about Bernstein and complete Bernstein functions and their relations to 
completely monotone functions. 

Definition 3.1 

A function / on R + is said to be completely monotone (CM) if it is infinitely 
differentiable and satishes the infinite set of inequalities: 

(—1)" D" f(y) > y > 0, for all non-negative integer n 

It follows from the definition and the Leibniz formula that the product of two 
CM functions is CM. A CM function can have a singularity at 0. 

Definition 3.2 

A function f on M. + is said to be locally integrable completely monotone (LICM) 
if it is CM and integrable over the segment ]0, 1]. 

Every CM function / is the Laplace transform of a positive Radon measure: 

Theorem 3.2 (Bernstein's theorem, (author?) j211) 

/(*)=/ e-'-V(dr) (9) 

^[0,oo[ 
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It is easy to show that / is a LICM if the Radon measure fi satisfies the inequality 



Theorem 3.3 

(i) If a,b > and the functions /, g are CM then the function a f + bg is CM. 

(ii) If f fa k = 1,2, ... is & sequence of CM functions pointwise converging to 



The proof of the second statement is based on Bernstein's theorem (Theo- 
rem [32]). Using Theorem 13. 31 a very useful theorem can be proved: 

Theorem 3.4 

If 

1. m is a positive measure on an interval I, 

2. the functions f(-,r) are defined and CM for all r in an interval I except 
perhaps for a set of zero measure m, and 

3. the integral g(x) :— J. f(x, r) m(dr) exists for x > 0, 
then the function g(x), x > 0, is CM. 

The last theorem can be proved by approximating the integral by finite Riemann 
sums and applying Theorem 13.31 

Definition 3.3 

A function f on R + is said to be a Bernstein function (BF) if it is non-negative, 
differentiable and its derivative is a CM function. 

Since a BF is non-negative and non-decreasing, it has a finite limit at 0. It can 
therefore be extended to a function on R+. 

Theorem 3.5 fJ^,\M] 

If /, g are CM then the pointwise product f g is CM. 

Let / be a Bernstein function. Since the derivative D/ of / is LICM, Bern- 
stein's theorem can be applied. Upon integration the following integral repre- 
sentation of a general Bernstein function / is obtained: 



where a > 0, b = D/(0) > 0, and u(dr) := n(dr)/r is a positive Radon measure 
on R + satisfying the inequality 




(10) 



/ then f is CM0/. 




]0,oo[ 




(11) 




(12) 
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The constants a,b and the Radon measure v are uniquely determined by the 
function /. 

Theorem 3.6 % [H [Hy 

(i) If f is a CM function, g is a BF and g(y) > for y > then the function 
fog is CM. 

(ii) f is a BF if and only if for every A > the function exp(— A f{x)) is CM. 
Corollary 3.7 % [I| 0/ 

If g is a non-zero BF then 1/g is a CM function. 

Note that the converse statement is not true: the function f(y) := exp(— y) 
is CM but 1// is not a BF. 

Definition 3.4 

A function f is said to be a complete Bernstein function (CBF) if there is a 
Bernstein function g such that f(y) = y 2 g(y). 

Theorem 3.8 fid] 

A function f is a CBF if and only if it satisfies the following two conditions: 

(i) f admits an analytic continuation f(z) to the upper complex half-plane; 
f(z) is holomorphic and satisfies the inequality Im/(z) > for Imz > 0; 



The derivative Dg of the Bernstein function g is a LICM function h. Hence 
we have the following theorem: 

Theorem 3.9 YlK] Every CBF f can be expressed in the form 



where h is LICM and a — g(0) > 0. Conversely, for every LICM function h 
and a > the function f given by (|13[) is a CBF. 

Proof. For the first part, let g be the BF in Definition 13.41 and let h := Dg. 
Since f Q h(x) dx = g(l) — g(0) < oo, the function h is LICM. For the second part, 
note that if h is LICM, then g(y) = a + h(s) ds is a BF and f(y) — y 2 g(y). 



(n) f{y) >0fory&R+. 



f(y) = vKv) + a v 



(13) 



□ 



Since the Laplace transform of a LICM function h has the form 




(14) 
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where fi is the Radon measure associated with h, every CBF function / has the 
following integral representation 

f(y) = b + ay + y[ ^ (15) 
7]0,oo[ r + v 

with arbitrary a,b = /i({0}) > and an arbitrary positive Radon measure /x 
satisfying eq. (|10[) . The constants a,b and the Radon measure /i are uniquely 
determined by the function /. 

Noting that y/(y + r) = r[l/r — l/(y + r)], we can also express the CBF / 
in the following form 

f(y) = b + ay + [l - e" z »] h(z) dz 
Jo 

where h(z) :— Jj , e~ rz m(dr) > and m(dr) := r /U(dr) satisfies the inequal- 
ity 

m(dr) 



[0,oo[ r ( r + 1) 



< OO 



Let v(dz) := h(z) dz. We have 

zv(dz) f°°zh(z)dz 



[o,oo[ ! + z Jo l + z J[o 



m(dr) [1/r - eT(0,r)] 



Using the asymptotic properties of the incomplete Gamma function [l[ it is 
possible to prove that the right-hand side is finite, hence the Radon measure 
is(dz) :— h(z)dz satisfies inequality (|12p. We have thus proved an important 
theorem: 

Theorem 3.10 

Every CBF is a BF. 

However 1 — exp(— y) is a BF but not a CBF. 

The sets of LICM functions, BFs and CBFs will be denoted by £, 03 and 
respectively. 

The simplest example of a CBF is 

y* oo 

Mv) ■= y/(y + *) = y 2 ^ sy [i - ^ sa ] ^ 

Jo 

a > 0. It follows from eq. (fT5|) that every CBF / which satisfies the conditions 
/(0) = and lim y ^oo f(y)/y = is an integral superposition of the functions 
(f a . The CBF ip a corresponds to a Debye element defined by the relaxation 
function G a (t) = exp(— at). 

We shall need the following properties of CBFs: 

Theorem 3.11 £H, [H/ 



7 



(i) f is a CBF if and only if yj f(y) is a CBF; 

(ii) if f,g are CBFs, then f o g is a CBF. 

The second statement follows easily from Theorem l3.8l 

Theorem 3.12 

y a is a CBF if < a < 1. 

Proof. 



Using Theorem 13.1 If ii). the last theorem implies the following corollary: 
Corollary 3.13 

Iff is a CBF and < a < 1, then f{x) a is a CBF. 

Remark. 3.14 /// is CM then f n is CM for positive integer n f Theorem ] 3. 5)) . 
However f a need not be CM for non-integer a. 

If f G 5 then f a need not be a CBF for a > (including integer values) unless 
< a < 1. 

4 Positivity of solutions in one- and three-dimensional 
space. 

Applying the results of the previous section, we get the following result: 
Theorem 4.1 

If a > and the relaxation modulus G is CM then the function Q defined by 
eq. © is a CBF. 

The mapping (a, G) £ R + x5-»Qe£ defined by eq. §5§ is bijective. 




and thus 




□ 



8 



A one-dimensional solution of eq.(U) is given by 

= U {1 \p, \x\) := A(p) exp(-B(p) \x\) 

with B(p) = fV 2 p/Q(p) 1 ' 2 and = ff (p)/[2B(p)]. If Q € ff, then Q(y) 1 / 2 
is a composition of two CBFs, namely y 1 ^ 2 (Theorem 13. 12[) and Q, hence it is 
a CBF by Theorem [3J3J The function B(p) is a CBF by Theorem EHU and 



1/B(p) is a CM function by Theorem 13 . 1 01 and Corollary 13. 71 

The amplitude of the solution of Problem I is given by A(p) = l/[2 pB(p)] + 

ua/[2B(p)]. The first term is a CM function because it is the product of two 

CM functions. The second term is also CM, hence A(p) is CM. The amplitude 

of the solution of Problem II A(p) = l/[2 5(p)] is also CM. 

For every fixed x the function exp(— B(p) \x\) is the composition of a CBF 

and the function B, which is a CBF and therefore a BF. By Theorem 13.61 the 

function exp(— B(-) \x\) is CM. This proves that for d = 1 the solutions of 

Problem II and Problem I with uq = are non-negative. 

In a three-dimensional space the solution u^(p, x) = U^'ijp, \x\) of (|4]) is 

given by the equation 



15j so that 



r,(3) 



1 



(p, x) = —A(p) B(p) eM-B( P ) \x\) (16) 
Anr 

But A(p) B(p) = g(p)/2. If u = then g is CM. Hence u^(-, x) is the product 
of two CM functions and thus CM. 



5 Positivity of scalar Green's functions in arbi- 
trary dimension. 

In an arbitrary dimension d 

^-j&mf^ MM+w ** <17) 

The above formula can be expressed in terms of MacDonald functions by using 
eq. (3) in Sec. 3.2.8 of (author?) pjj|: 

-d/4-1/2 „C_\ n d/2-l 

where B(p) is defined in the preceding section. 

The MacDonald function is given by the integral representation 

poo 

Kp(z) = / exp(— z cosh(s)) cosh(/^s) ds (19) 
Jo 
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Since cosh(?/) is a positive increasing function, it follows immediately that K^z) 
is a CM function. 

We shall need a stronger theorem on complete monotonicity of MacDonald 
functions. 

Theorem 5.1 fM]- 
The function L^(z) ■= z 1/2 K^z) is CM for n > 1/2. Q 

The proof of this theorem requires a lemma. 
Lemma 5.2 

If a > then the function (1 + l/x) a is CM. 

Proof. We begin with < a < 1. Setting t = l/(xy) we have that 

o_ r°° dy f 1/x t*- 1 _ 

x a Ji y 1+a (xy + iy- a a J (1/i+l) 1 -" 

rl+l/x f ■ 

a J u a - 1 du= I 1 + - 

Since for each fixed value of y > the function (xy + l)"^ 1 is CM, the function 
(1 + l/x) a (x > 0) is also CM. 

The function 1 + l/x is CM, hence for every positive integer n the function 
(1 + l/x)" is CM. We can now decompose any positive non-integer a into the 
sum a = n + (3, where n is a positive integer and < /3 < 1. Consequently 

(1 + l/x) a = (1 + l/x)" (1 + l/a;)^ 

is CM because it is a product of two CM functions. 

□ 

Proof of the theorem. For /i > —1/2 the function L^(z) has the following 
integral representation: 

(20) 

11] , 8:432:8. By Lemma [5.21 the integrand of the integral on the right-hand 



side is CM if /i > 1/2. Hence the integral is the limit of sums of CM functions, 
therefore itself a CM function. Consequently, the function L^(z) is the product 
of two CM functions. By Theorem [33] is CM. □ 

We now note that u {d \p, x) = C p( d ~ 3 )/ 2 g(p) F[p), where C = p^- 1 / 2 / \2n) d / 2 
and F{p) = Q{v)~ {d+l)/i Ld/2-i{B(p) r). We shall prove that F(p) is the prod- 
uct of two CM functions. 



1 The theorem is valid for fj, > 0, see 12111 . but we do not need this fact. 
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Lemma 5.3 

IfQ is a CBF and a > then Q(p)~ Q is CM. 

Proof. Let n be the integer part of a, a = n + f3, > j3 < 1. is 
CM (by Theorem [330] and Corollary ELZ} and therefore also <3(p)~™ is CM. By 
Theorem 13 . 1 1 1 the function Q^p) 13 is a CBF, hence l/Q(p)P is CM. Consequently 
Q(p)- Q is CM. 

□ 



Lemma 5.4 

IfQeg and /i > 1/2 then L d/2 _ 1 (B{p) r) is CM. 

Proof. We have already proved that B is a BF. It follows from Theorem 15.11 
and Theorem 13.61 that the function L fl (B(p) r) is CM for \i > 1/2. 

□ 

Lemma IS~3l implies that for every d > the factor Q(p)~( d+1 ^ 4 is CM while 
Lemma I5T41 implies that L d / 2 -\{B(p) r) is a CM function of p for every r > for 
every d > 1. Hence the function F(p) is CM. The factor p( d - 3 )/ 2 is CM if d < 3. 
Consequently, for d < 3 the solution a;) of Problem II is non-negative. The 
solution of the same problem with an arbitrary source of the form s(t) 5(x) and 
s(t) > can be obtained by a convolution of two non- negative functions and 
therefore is also non-negative. 

For d < 5 Problem I with uq = uq = has a non-negative solution if c\ > 0. 
For d < 7 Problem I has a non-negative solution if c\ = and c 2 > 0. 

We summarize these results in a theorem. 

Theorem 5.5 

In a viscoelastic medium of dimension d < 3 with a constitutive relation 

a = ae + G(t)*e, a > 0; Ge5 

Problem II as well as Problem I with the initial condition uq = have non- 
negative solutions. 

For zero initial data Problem I has a non-negative solution if d < 5 and 
ci > 0, or if d < 7, c\ = and c 2 > 0. 

Anisotropic effects can be introduced by replacing the operator V 2 by g kl dk di . 
If h H g lm = 8 k m then 



n d/i-l/2 n ( n \ n d/2-l 



,«*/4-l/2 ^(p) pd/2-1 
(27t) d / 2 Q(p) d / 4 + 1 /2'' 

itya-i (p 1/2 WQ(p) 1/2 ) (21) 



If Q is a CBF then u(t, x) > 0, with r := [h u x k x l ] 1/2 , cf (author?) 
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6 Positivity properties of vector-valued fields. 



It is interesting to examine the implications of CM relaxation kernels on positiv- 
ity properties of vector fields. We shall prove that in a simple model complete 
monotonicity of a relaxation kernel implies that the Green's function is positive 
semi-definite. 

Unfortunately the tools developed in Sec. [3] fail for matrix- valued CM and 
complete Bernstein functions. In particular, the product of two non-commuting 
matrix- valued functions need not be a CM function and the function / o G, 
where / is CM and G is a matrix-valued BF, need not be CM. 

Definition 6.1 

A matrix-valued function F : M + — > R nxn is said to be a CM function if it is 
infinitely differentiable and the matrices (—1)™ D™ F(y) are positive semi-definite 
for all y > 0. 

Definition 6.2 

A matrix-valued Radon measure M is said to be positive if 



for every vector v 6 M™ and every non-negative function f on R+ with compact 
support. 

It is convenient to eliminate matrix-valued Radon measures by applying the 
following lemma [HI ]: 

Lemma 6.1 

Every matrix-valued Radon measure M has the form M(dx) = K(x)m(dx), 
where m is a positive Radon measure, while K is a matrix-valued function 
defined, bounded and positive semi-definite on R + except on a subset E such 
that m(E) = 0. 

Theorem 6.2 fU l 

A matrix-valued function F : M.+ — * R™ x ™ is CM if and only if it is the Laplace 
transform of a positive matrix-valued Radon measure. 

The following corollary will be applied to Green's functions: 
Corollary 6.3 

If R(p) := f e~ pt R(f ) dt is a matrix-valued CM function then R(i) is positive 
semi-definite for t > 0. 

Definition 6.3 

A matrix-valued function G : R + — » Ml™*™ is said to be a Bernstein function 
(BF) if G(y) is differentiable and positive semi-definite for all y > and its 
derivative DG is CM. 
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Definition 6.4 

A matrix-valued function H : M. + — > R™ xn is said to be a complete Bernstein 
function (CBF) if H(y) = y 2 G(y), where G is an n x n matrix-valued BF. 



The integral representation (|15p of a CBF remains valid except that the 
Radon measure has to be replaced by a positive matrix-valued Radon measure 
N(dr) = K(r) i/(dr): 



H(y) = B + yA + y / 

•/10,c 



i )0O [ r + y 

where the Radon measure ^ satisfies the inequality 

^(dr 



K(r) i/(dr) , N 

(22) 



[0,oo[ !+ r 



< oo (23) 



the matrix- valued function K(r) is positive semi-definite and bounded z^- almost 
everywhere on R+ while A, B are two positive semi-definite matrices. Every 
matrix-valued CBF H can be expressed in the form 

H(y)=yF(y) + yA (24) 

where F is a matrix-valued LICM function. 

We now consider the following problem for the Green's function Q 

pD 2 g = ADQ + G *V 2 DQ + 8(t)5(x)I, t>0, x e R d (25) 

where A is a positive semi-definite nxn matrix and G is an n x n matrix- valued 
relaxation modulus. 

If the relaxation modulus G is a CM matrix-valued function then the func- 
tion Q(p) := pG(p) is a matrix-valued CBF. The function Q is real and positive 
semi-definite, hence it is symmetric and has n eigenvalues qi (p) and n eigenvec- 
tors &i, i = 1, . . . , n. We shall now assume that the eigenvectors are constant: 

n 

Q(p) = ^2<li(p)ei ® e * 

i=l 

It is easy to see that the functions qi, i — 1, . . . , n, are CBFs. 

The Laplace transform Q{p, x) of the Green's function is given by the formula 

G(p, x ) = T^fj ei(Kx) [p 2 1 + l k l 2 Q(p)] 1 d * k = 

£ j^- d J ^' X) [P 2 + |k| 2 <fc(p)] e, ® ei d d k ee 

d 

Qi(p,x) et ® Bi 



i=l 
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where 

d/4-1/2 d/2-1 

and Bj(p) = p 1 / 2 p/qiip) 1 ^ 2 , i = 1, ...,n. Assume for definiteness that d < 
3. The argument of Sec. [5] now leads to the conclusion that the functions gi, 
i= 1, . . . , n, are CM. Consequently the Laplace transform Q(-, x) of the Green's 
function is a matrix-valued CM function for x G M, d and the Green's function 
Q(t,x) is positive semi-definite for t > 0, x G K d . In particular, we have the 
following theorem: 

Theorem 6.4 

Let p G R+, d < 3, s(t, x) = S(t) 8(x) w, where w G 1". 

If G(s) = X)"=i Gj(s) ej (g e^ and A = a « e » ® e « with CM functions 

Gk and reai numbers > 0, i = 1, . . . , n, then the solution u of the problem 

pD 2 u = ADu + G*V 2 Du + s(t, x), t>0, x G M d 

satisfies the inequality 

(u(t, x), w) > 0, < > 0, x G M d . (26) 

7 Positivity in isotropic viscoelasticity. 

Consider now the Green's function Q of a 3D isotropic viscoelastic medium. The 
function Q is the solution of the initial-value problem: 

pD 2 Gkr(t, x) = G(t)klmn * D Q mr , n i + S(t) S(x) 5 kri 

i>0, i£l 3 , fe,r = l,2,3 (27) 

with zero initial conditions, and 

Gklmn{t) = A(t) (28) 

where the kernels X(t),p,(t) are CM and p G R+. The function G with the 
components Gkimn takes values in the linear space & of symmetric operators 
on the space S of symmetric 3x3 matrices. It is easy to see that under our 
hypotheses this function is CM: 

(-1)" (ei, G(t) e 2 ) > for all n = 0, 1, 2 . . . 

for every ei, e2 G S, where (v, w) := v k i w k i is the inner product on S. 
The Laplace transform Q of Q is given by the formula 

g( P ,x) = — {v®va- 1 f l (p,|x|)+ 
pp L 

[I-V0VA- 1 ] F T (p,\x\)} 
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where A := V 2 , 

\2 



Fh{pjr) ^ S jjPL c -p 1/2 s,( P )r (2g) 

F^p^-S-lM^e-P^sArtr (30) 



and 



tW :- JL (32) 

Since <zl(p) = p/sl(p) andqT(p) =p/st(p) are CBFs, the functions p/qi^ij)) 1 ^ 2 = 
p 1 ^ 2 sl(p) and p/qrip) 1 ^ 2 — p 1 ^ 2 st(p) are BFs. Hence the exponentials in 
eqs ([29l l30| are CM functions of p. Moreover the functions p sl(p)~ 2 , p st(p) -2 
are CBFs and thus sl(p) 2 /p and st(p) 2 /p are CM. It follows that the functions 
Fl(p, f) and Fr(p, r) are CM and therefore they are Laplace transforms of non- 
negative functions F^(t,r) and Fx^r). Their indefinite integrals ft,(t,r) := 
Jo r) ds and /t(£, r) := J -Fl(s, r) ds are also non-negative. The functions 
hi,(t, r) := A -1 fh(t, r) ds, /it(£j r) = A -1 /t(£, r) involve a convolution with a 
non-negative kernel and therefore are non-negative. The Green's function can 
be expressed in terms of these functions: 

G{t,x) = ~{V(g>Vfc L (t,|a:|) + [AI-V®V] h T (t,\x\)} (33) 

We shall use the notation v > if Vk > for k = 1,2,3. 
Theorem 7.1 

Let u = Vtfi + V x ip be the solution of the initial-value problem 

pD 2 u = G * V 2 D u + s(t, x), t>0, xeU d (34) 

with u(0, x) = = Du(0, x) and s(t, x) = Vf(t, x) + V x g(i, x). 

Then Vf{t,x) > for ail < > 0, x 6 M 3 implies that V(f>{t,x) > for all 
(>0,i£ M 3 . 

Similarly, V x g(t, x) > for all ( > 0, i £ I 3 implies that V x x) > 
for all i > 0, x £ M 3 . 

Proof. Substitute u = V</> + V x ip, s = V/ + V x g in the formula 

u(t,x)=J J G(t - s,x - y)s(s,y)d 3 yds 

where Q is given by (|33p . Noting that A -1 is a convolution operator commuting 
with V and V A -1 div s = V/ we have 



V<f>(t,x) = - f h(t- 
P J 



s, |x - y\) (V/)(s,y) d 3 y 
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We now note that [I - A -1 V ® V] s = s-V/ = Vxg. Hence 

V x i/)(t, x) = - / f T (t -s,\x- y\) V x g(s, y) d 3 y 
P J 

The functions /l and /t are non-negative, hence the thesis follows. 

□ 



8 Generalized Cole-Cole relaxation modulus. 

A particular example of a CBF is the rational function F{p) = Rm(j>)/ Sm{p), 
where Rm and Sm are two polynomials with simple negative roots k = 
1, . . . ,N, fii, I — 1, . . . , M , M — N or N + 1 satisfying the intertwining condi- 
tions: 

< Ai < pi < . . . fi N [ < Ajv+i ] 

(the last inequality is applicable only if M = N + 1). The function F is thus a 
good example for the function Q. 

A more general CBF is obtained by substituting in F the CBF p a , with 
< a < 1: 

F a ( P ) =R N (p a )/S M (p a ) 
(Theorem 13. 11 1) . The choice of Q — F a corresponds to a generalized Cole-Cole 
model of relaxation. For N = M = 1 the original Cole-Cole model 0, 13, EE E3] 
is recovered. 

Another useful CBF is the function F T , a (p) := F ((1 + (rp)) Q ), for some 
a e [0, 1] and r > 0. The function (1 + {rp)) a is a CBF, hence F T a G J too. 
This function has a low-frequency behavior appropriate for poroacoustic and 
poroelastic media. A simple poro-acoustic model of this type is defined by the 
equation^, M, S S 0, M, Q i 

pp 2 K(p) u(p, x) = A t2(p, x) + s(p, x) (35) 

where K (p) := F T , a (p)~ 1 is a CM function, known as the viscodynamic operator[2: 
Bernstein's theorem implies that the inverse Laplace transform H (t) of K(p) is 
non-negative. In terms of the kernel H(t) eq. (133)) is equivalent to the time- 
domain equation 

p D 2 H(t) * u(t, x) = A ie) + s(i, x) (36) 

The results of Sec. [5] imply that for Q = F, F a , F a ^ T , and < a < 1, r > 0, 
d < 3 the Green's functions of {T]) and (f36|) are non-negative for d < 3. 

9 Beyond viscoelasticity. 

9.1 Positivity for Q 1 ' 2 £ J. 

The proofs of positivity in the preceding sections were motivated by the con- 
stitutive equations of viscoelasticity which led to the assumption Q s 5^ and its 
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variants underlying all the proofs. We now approach the same problem from 
the other end and ask whether positivity would hold under weaker assumptions 
than Q G 

Positivity can be shown to hold under somewhat weaker assumptions. The 
Green's function for eq. |T]) can be expressed in the form 



{2n) d / 2 p 



(d-3)/2 



B(P) 
P 



(d+l)/2 



r -(d-i)/2 L {B {p)r) (37) 



For d = 1,3 the function £d/2-i( z ) = (tc/2) 1 ^ 2 exp(— z) and the last factor in 
(|3"T|) is CM for every r if and only if B is a BF (Theorem I3.Q|) . On the other 
hand, if B is a BF then, by Lemma [A. 11 the function B(p)/p is CM. If d is 
an odd integer then the third factor on the right-hand side of (f3"T|) is CM. The 
second factor is CM if d < 3. Therefore the Green's function is non- negative if 
B is a BF and d = 1,3. 

Positivity of can be proved under somewhat stronger assumption B G 5. 
Let < 7 < 1. Since -B,p 7 € $, from Theorem IB.2I we conclude that 



P 



B(p) 



1/2 



n 7/2 



Hence 



Mp) 



B(p) 



P 



3/2 



,-7/2 _ 



B(p) 



P 



B(p) 



P 



1/2 



'7/2 



is the product of two CM functions. Hence it is CM. The pointwise limit 
lim 7 ^i/ 7 (p) is CM. Since L (B(p)r) is CM, eq. (|37|) implies that U^(-,r) 
is CM for all r > 0. Hence u^(t, x) > 0. 

If B £ # then the function Q{p) 1 ^ 2 = p/B(p) is also a CBF but Q need not 
be a CBF. Since Q G 5 entails Q 1 ^ 2 G 3 but not conversely, the assumption 
that B G 5 or even B G 05 is weaker than the assumption Q G 3- For d = 1 
and 3 this constitutes a generalization of the results obtained in Sec. O 



9.2 Example. 

We shall consider viscoelastic media with power law attenuation, defined by the 
formula 0, EH 

Q(p)= M , ° a _ lia , a>0, c >0, 0<a<l (38) 
[1 + ap Q 1 J Z 

Time-domain relaxation calculus and creep compliance for power law attenu- 
ation can be found in (author?) [IBj]. Applications of this model in acoustic 
inversion and seismic inversion can be found in (author?) [25| and (author?) 
(lij . For d = 1,3 the Green's function can be expressed in terms of an expo- 
nential 

U {d \p,r) = F^ d \p,r)e-^ r) 
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where F^' (p, r) is an algebraic function and the phase function 

$(p,r) =p 1/2 pr/Q( P y/ 2 =p 1/2 [p + ap a ] r/c (39) 
The attenuation is non-negative: 

| e -»(-iu,,r)| ^— R.c — icj,r) g—ar a uj a cos(7ta/2) ^ 

and the Green's function is causal by a Paley- Wiener theorem ((author?) [23j ]. 
Theorem XII) because 

|ln|F( d )(r,p)e- <I, (- iw ' r )|| , a p x l 2 r cos(W2) f°° oj a , 

1 ■ dw w — — / duj < oo 



1 + w 2 c jj 1 + 

The above results can be extended to d — 2 by using the asymptotic formula 
K (z) ~ ^n/2z exp(-z). 

Theorem 9.1 

(i) The function Q 1 / 2 is a CBF for < a < 1. 

(ii) The function Q is a CBF for 1/2 < a < 1. 

Proof. We note that for positive real p both Q(p) and Q(p) 1 ^ 2 are real and 
non-negative. If < argp < 7t then 

* := arg [l + ap 11 " 1 ] > arg [ap a_1 ] = -(1 - a) argp > -(1 - a) n 

and * < 0. Thus 

< arg [1 + ap a - l \ ~ X < (1 - a) n 
By Theorem GHQ 1 / 2 S However 

< arg [l + ap^ 1 ]^ < 2(1 - a)n < 71 
holds only for a > 1/2. Hence Q is a CBF for 1/2 < a < 1. 

□ 

The results of the previous subsection imply that the Green's function is non- 
negative for < a < 1. 

Corollary 9.2 

The relaxation modulus G is CM if and only if 1/2 < a < 1. 

Proof. G is CM if and only if Q is a CBF. Hence G is CM for 1/2 < a < 1. 

We shall prove by direct calculus that G is not CM and therefore Q is not a 
CBF for < a < 1/2. The relaxation modulus is given by the inverse Laplace 
transform 

G(t) = —[e pt — —dp (40) 

V ; 2niJ ts pil + ap"- 1 } 2 y y ' 
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Figure 1: Contour deformation. 



where B denotes the Bromwich contour running parallel to the imaginary axis 
in the right-half complex plane. For a > and < a < 1 the denominator does 
not have zeros on the Riemann sheet —tc < argp < tc and all the singularities of 
the integrand are situated on the cut along the negative real axis. The integrand 
does not have a pole at nor on the cut. For t > the Bromwich contour can 
be replaced by a half-circle at infinity in the half- plane Rep < and the Hankel 
loop H encircling the cut from — oo below the cut to zero and back to -co above 
the cut (Fig. [I]). The contribution of the half-circle vanishes, hence 

c 2 f 1 

G(t) = — / e pt ? dp = 

2tti J h p [1 + apa-i] 2 

c 2 f°° 1 

— / e~ flt Im — F ; — r . — dR = 

Wo R[l + a(R exp^TT))"- 1 ] 2 



c 2 f°° _ Rt 2a R a sin(Tta) - a 1 R 2a ^ sm(n(l - 2a)) 



7t Jo [R + a 2 R 2 *- 1 - 2a R a cos(TTa)] 2 



dR 



The integrand is non-negative for a > 1/2 and changes sign for a < 1/2. By 
Bernstein's theorem the relaxation modulus G is CM for a > 1/2 and is not 
CM for a < 1/2. Fig. [2] shows that G is neither monotone nor non-negative for 
a < 1/2. 

□ 

It follows that Q is not a CBF if a < 1/2. Hence we have extended positivity 
to the case when Q g $ but Q 1 / 2 e 



19 







5 







4 







3 







2 







1 




20 



40 



60 



80 



100 







1 



Figure 2: The relaxation modulus for the power-law attenuation. The exponent 
values are a = 0.3, 0.4, 0.5, 0.6, 0.9, from bottom to top. 

10 Concluding remarks. 

In an elastic medium and in the case of Newtonian viscosity a non-negative 
source term excites a non-negative displacement pulse. This positivity property 
is shared by hereditary viscoelastic media with non-negative relaxation spectra 
and is limited to 1,2 and 3-dimensional problems. The concept of positivity 
can be extended to systems of partial differential equations (vector fields). For 
general matrix-valued pseudo-differential operators (anisotropic viscoelasticity) 
positivity holds in a non-local sense. In isotropic viscoelasticity it holds for 
longitudinal and transverse waves separately. 

The example in Sec. 19.21 shows that positive monotonicity of G is not a 
necessary condition for positivity. In Sec. 19.11 positivity is shown to hold for 
partial differential equations of a more general type for dimension < 3. 
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A A lemma. 

Lemma A.l 

If f is a BF, then the function f(x)/x is CM. 
Proof. Let 



J]0,oo[ 

with a, b > and a positive Radon measure m. By Bernstein's theorem the 
function 



is the superposition of CM functions with non-negative weights. By Theorem l3.4l 
the function f(x)/x is CM. 






is CM. Therefore 




□ 
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B Some properties of Stieltjes functions and CBFs. 

Definition B.l 

A real function f on R + is a Stieltjes function if it has the integral representation 



with a > and a positive Radon measure [i satisfying (|10p . 

The Stieltjes functions form a cone 6: 
if a, b > and f,g e 6 then a f + bg e &. 
Every Stieltjes function is a CM function. 

Theorem B.l 

A function f is a CBF if and only if f(x)/x is a Stieltjes function. 

Theorem IB . 1 1 follows from the integral representation (fT5|) . Compare Theo- 
rem [Bj] with Corollary 13. 71 

We shall also need the following theorem: 

Theorem B.2 

If f, g G $ and < a < 1 then h := f a g x - a e $. 

Proof. If /, g G 5 then / and g map R+ into R + and C+ into itself. Conse- 
quently h maps R+ into IR+ and 




< arg h(z) <a7t+(l — a)n = n 



for z G C+. Hence, by Theorem I3.8[ ft e 5- 



□ 
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